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'. 1. Introduction and statement of the theorem 
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The simplest model of scattering/quantum resonances comes from considering compactly 
supported potentials on the real line, 

(1) V{x) e M, \V{x)\ < C, V{x) = for |x| > L , 

and the corresponding Schrodinger operators, 



m 
(N 

> 
m 

O : (2) Hy = -dl + V{x) 

(N ■ 

O ' on M, or on [0, oo) with Dirichlet or Neumann boundary conditions. 

' 

Q \ The resonances or scattering poles of Hy are defined as the poles of the meromorphic 

^ ' continuation of the resolvent, Rv{^) = {Hy — from ImA > 0, to C. Except for 

CLc the poles at A for which A^ are eigenvalues of Hy, -Ry(A) is bounded on for ImA > 
fH ! 0. Its Schwartz kernel, that is the Green function, continues meromorphically across the 
continuous spectrum corresponding to M. Its poles are the resonances of Hy. 

An illustration based on the numerical codes of [1] is given in Fig{T] The poles on the 
positive imaginary axis correspond to the bound states oi Hy, and the poles on the negative 
are called antibound states. Note that they appear to be exactly symmetric with the bound 
states. In this note we prove a simple theorem inspired by numerical experiments using [1]: 



Theorem. Consider the Dirichlet (or Neumann) boundary condition on [0, oo) and a com- 
pactly supported piecewise continuous potential Vq, supp Vq C [0, A). Let Vi > 0, B > A, 
and put 

Vix)=Vo{x) + l[A,B]i^Wl- 

Then the bound and antibound states of H^iy with moduli greater than some fixed /cq > 
are symmetric modulo errors of size e~'^'^, q ^ oo. 

Equivalently we can consider the semiclassical problem 

i-{hd,y + Vix))u{x) = z{h)u{x) , 

for which the conclusion of the theorem says that bound and antibound states with moduli 
greater than hk^ are symmetric modulo exponentially small errors, exp(— c//i), as h 0. 
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Figure 1. The (color coded) correspondence between classical dynamics 
and the distribution of resonances. The potential and its resonances are 
obtained using splinepot (40* [0, 1 , -2 , 1 , 0] , [-2 , -1 ,0 , 1 ,2] ) from [1]. 



We think of as a barrier separating the potential in the interaction region, 

Vo(x), from infinity. The same results hold on the line but the proof becomes slightly 
more cumbersome to write. In Fig|2]we show an example of potentials Vq, W = Vi1[a,b], 
and V = Vq + W. It is quite possible that using more sophisticated methods - see for 
instance [10] and [15] - more general barriers can be considered. Our goal here was to 
present a simple new result discovered by a numerical observation. It is easy to see ( [T7] 
and ^ below) that for the problem on the half line the bound and antibound states are 
never exactly symmetric. Yet, in a presence of a mild barrier, they are symmetric within 
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numerical accuracy of a computation: the exponential convergence is indeed very rapid. 
This is illustrated in FigjS) we plots of imaginary parts of bounds states and negatives 
of the imaginary parts of antibound states for q'^V{x) as a function of g^. The difference 
between the two pictures is striking. As expected the bound states are not much affected 
by the presence of W but the antibound states change dramatically and in the presence of 
a barrier become nearly symmetric (this is a curious pseudospectral effect). The high lying 
states of Vo also exhibit the symmetry - experiments show that it is always so, even when 
there is no barrier, and it improves for more regular potentials. 

The study of resonances/scattering poles in one dimension has a long tradition going 
back to origins of quantum mechanics, see for instance [13] . Perhaps the first study of their 
distribution was conducted by Regge [16]. For mathematical results in one dimension see 
[I] , [S| , [S| , [H] , [El , [IZ] , [20] , and many other articles. Concerning antibound states, Hitrik [H] 
showed (using a Ricatti equation approach which we also find useful in §2]) that for positive 
compactly supported potentials, there are no antibound states in the semiclassical limit. 
That of course corresponds to our result since there are no bound states either. Simon 
[T7] showed that for a half line problem existence of n bound states implies the existence 
of n — 1 antibound states. Since the set of resonances of an even potential is the union 
of Dirichlet and Neumann resonances of the half line problem, this means that having n 
bounds states implies the existence oi 2n — 2 antibound states. As can be checked using 
[1] this is often optimal for negative potentials but never for potentials with a barrier. 

Our note is organized as follows: in §2] we give the elementary proof of the theorem and 
in ^ we describe the ideas behind the computation of resonances in one dimensions. The 
MATLAB codes based on that section are available at [1]. 

Acknowledgments. The work of the second author was supported in part by a National 
Science Foundation grant DMS-0200732. 

2. Proof of the theorem 

We will prove the theorem for A = 1 and B = 2 - the general case is identical. We 
consider transfer operators for solutions of Hg2Y + k"^: 

Mo{k):[u{0),u'{0)]^[u{l),u'{l)] 

which depends only on Vq, q, and k, and 

Mi(A;):[w(l),w'(l)]^[w(2),w'(2)]], 

which is completely explicit since we know Vi. 

^^'f*) - ^ ( |:-„t ^"h*. ) • > ^ + '/^- 

As in the code described in ^ ik, k > 0, is a bound state if and only if 

Mi{k) o Mo{k) [0, 1] = [a, -ka] , 
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Figure 2. Potentials with bounds and antibound states shown in Figj3l 



for some a, and —ik is an antibound state if and only if 

M,ik)oMoik) [0,l] = [(3,k(3], 

for some /5 (note that same k will never do for both, so they are nei'er exactly symmetric). 

The conditions for a bound (— ) and an antibound state (+) then become (note that the 
left hand side depends on k and q only, and not on ±): 

^ ^ _ l-/j±(fc)exp(-2fci) ^ 1 ± k/k, 

^ ' u{l) 'l + /5±(A;)exp(-2A;i) ' ^ l^k/ki' 
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Figure 3. Imaginary parts, A, of bounds states (corresponding to eigenval- 
ues at — A^), and negatives of the imaginary parts of antibound states of g^Vo, 
and of q^{Vo + W), as functions of q"^. 

In fact, putting v{k) =^ we use ([3]) to obtain the following equations 

bound/antibound states: 

±^{h{l + e-^'^) + v{k)il - e-^'^O) = hil - e-''^) + v{k)il + e-^"^) , 
h 

or 

ki (^±^{1 + e"^"^) - (1 - e-2fc^)^ = v{k) (^1 + e-2fci ^ ^(i _ g^^/cr 
from which (HI) follows. 
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The behaviour of v{k) as k varies is well known: it is monotonic between — oo and 
oo where oo correspond to k"^ which are Dirichlet eigenvalues of Hg2Y^ on [0, 1]. But the 
equations for ± are the same up to exponentially small errors! 

More precisely, suppose that 

u" = {k^ + V{x))u, u{0) = 0, u{0) = l, 

u = u{x, k). Consider 

v{x, k) =^ k)/u{x, k) . 
Then v satisfies the Ricatti equation 

v' = e + V -v'^ . 

Denote differentiation with respect to k by i). We get the following equation for v: 

v' = 2k- 2vv, 

which we can solve by the method of integrating factor. Noting that u{x)'^v{x)\r^=Q = 
(from the boundary conditions) we get 

2k 

v{x,k) = — — — / u{y,kfdy 
u{x, ky Jo 

and in particular we get an expression for the derivative of the Dirichlet to Neumann map: 

2k 

ky Jo 

which we can estimate from below as follows. 

Since we assumed that supp V C [0, 1), for some e > we have 

u{x, k) = Ae''' + 5e-^^ , 1 - e < x < 1 

Then 

1 \ , a — 1 dcf ^ 2A- 

a + 1 B 

and using the same notation, 

2k 

v{^,k)>—-—— I u{x,kfdx 
u[l,ky J^_, 

(5) " + 

(a + 1)2 
>l-5, k>ko{e,5), 

for any 5 > 0. 
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We recall that the condition (j4]) for being a bound (— ) or an antibound (+) state was 

v{l, k±) = {kl + gVi)^(l + g±{k±, q)) 

where g± = (9(e~^'^). Put 

Fjk)'^ "(^-^^ , -1, 
{k^ + qWi)2 

so that F{k±) = g±{k±,q) = 0{e^'^'^). Once we show that F{k) ^ 0, we will know that 
the roots of F are stable, and by standard theory, small perturbations to the equation lead 
only to small perturbations to the roots. 

More precisely, we use (ED to estimate 

v{l,k){e + q^Vi) - kv{l,k) 



> 



(k^ + q^Vi)^ 
il-6)ia + lfie + q^Vi)-eia^ 
(P + g2y,)l(a + l)2 



>-, k>ko{e,6), 

provided that 6 is taken small enough depending on Vi > 0. Hence, by the mean value 
theorem there exists some < s < 1, such that 

\k -kJ- 



\F{{l-s)k+ + sk.)\ 
<Ce~^'''{{l-s)k+ + sk^)/S 
< e'"^ , k±> ko, q> qo. 

Note that we used the fact that k± are necessarily bounded by Cq. 

Replacing the explicit solutions by WKB approximations might give a more general 
result. 



3. Numerical computation of resonances in one dimension 



In this section we describe the ideas behind the codes, squarepot.m and splinepot.m, 
used to produce Fig{T]and[3l These MATLAB codes are available at |4]. 

If the support of V is contained in a compact interval [—L,L], we can compute both 
resonance solutions and ordinary eigenvalues of the Schrodinger problem, [Hy — \^)u = 0, 
by writing appropriate boundary conditions at ±L: 

(Hv - A2)n = for X G L), 
(6) {dx + iX)u = at a; = L, 

{dx — i\)u = at X = —L. 
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In terms of A, this is a quadratic eigenvalue problem. We can introduce a new variable 
ip = Xu to convert this problem to a linear eigenvalue problem in two fields: 



(7) 



HyU — Xip = 

Am - ^ = 
{dx - iX)u = 
{dx + iX)u = 








for X G {—L, L), 
for X G [—L, L], 

dbti X —— Lj 
dbti X —— L. 



We now discretize the boundary and domain operators to get a finite-dimensional general- 
ized eigenvalue problem. For small discretizations with up to a few hundred unknowns, we 
can solve this generalized eigenvalue problem using MATLAB's eig command, which uses 
the dense eigensolvers in LAPACK [2]. For larger discretizations, we use MATLAB's eigs 
to call ARPACK, a standard Arnoldi-based iterative eigensolver [T2] . 

For the calculations shown in this note, we used a high-order pseudospectral collocation 
method to discretize the operators [12], [S]. We partition the support interval [—L,L] into 
subintervals, and approximate m by a high-order polynomial on each subinterval. At the 
Chebyshev points on the interior of each subinterval, we insist that the domain differential 
equations be satisfied exactly, while at the junctions between neighboring intervals, we insist 
that the solution u and the first derivative dxU must both be continuous. Assuming that 
the potential is smooth except possibly at the endpoints of the subintervals, the collocation 
scheme we use is spectrally accurate; that is, the error asymptotically decreases faster than 
any algebraic function of the order of the collocation scheme. As a simple check on the 
accuracy of the computed eigenvalues of ([7]), we increase the order of the method by 50%, 
recompute the eigenvalues, and compare the results obtained from the coarser and the finer 
discretization. 

We can write the analogue of in higher dimensions, with a Dirichlet-to- Neumann 
(DtN) map - or some approximation to a DtN map - in place of the boundary conditions 
at ±L. In more than one space dimension, this boundary map ceases to be a linear function 
of A, and so we cannot easily convert the problem into a linear eigenvalue problem. Re- 
searchers are studying these more complicated nonlinear eigenvalue problems for a variety 
of engineering problems [3]. Many of these problems involve resonances in models of elastic, 
acoustic, or electromagnetic resonators with radiation losses. 

For comparison, we will also discuss other methods for computing resonances. They are 
essential for effective codes for higher dimensional problems for which analogues of ([6]) are 
unavailable or become more complicated. 

Often, resonances are computed by changing the equation so that it is no longer posed 
on all of M, but instead is posed on some interval (— M, M) with homogeneous Dirichlet 
or Neumann boundary conditions. For example, if the support of V lies strictly within 
the interval {—L,L), we might add a complex absorbing potential outside of {—L,L), or 
we might scale the coordinate system into the complex plane by the method of perfectly 
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matched layer^. The change to the equation should be designed so that the modified 
equation mimics the behaviour of the original problem in the range {—L, L). 

To be more concrete, suppose that we modify the equation on the interval (L, M) so 
that we still have a nonsingular, second-order, ordinary differential equation in x whose 
coefficients depend on A. Now we specify two linearly independent solutions 7_|_(a;, A) and 
7_(a;, A) on (L, M) which satisfy the modified domain equation together with the initial 
conditions 

/ON 7+(^, A) = 1, d^'j+iL, A) = iX 

^ ' 7-(^,A) = l, d^^.{L,X) = -tX. 

These initial conditions are consistent with the conditions for outgoing and incoming waves 
on (L — e, L). Now suppose that 7(x, A) satisfies the differential equation on (L, M), and 
also the boundary condition 7(M, A) = 0. Then 

(9) 7(a;,A) = c(7+(a;,A)+p7_(x,A)) 

where c is an arbitrary constant and 

def 7+(M,A) 
7_(M, A) 

is a constant whose amplitude reflects how well the equation on (L, M) serves to absorb 
outgoing waves. We can therefore convert the condition at x = M to a condition at x = L. 
Subsituting ([8]) into Qj, we have 



which, for regions of the complex plane where |p(A)| is small, can be treated as a pertur- 
bation of the exact outgoing wave condition at L. 

In summary, by changing the Schrodinger equation outside the interval {—L, L), imposing 
homogeneous Dirichlet boundary conditions at ±M, and then transporting the conditions 
at ±M to conditions at ±L, we arrive at the equations 

{Hv - A2)n = for X G L), 

(10) {dx + iX)u = at X = L, 

{dx — iX)u = at X = —L. 

where 



A = A 



P(A) \ 



,l + p(A)y 

For values of A where |p(A)| ^ 1, ([6]) and (fTOj) may be treated each as a perturbation 
of the other. We note that p(A) and /3-|-(/c) exp(— 2/ci) of (j4]) play similar roles in the two 
situations. However, the smallness of p(A) is achieved through ellipticity due to the complex 



"'^See [7] for a comparison of that method with the complex scahng method described, for instance, in 
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deformation, and the smallness of P±{k) exp(— 2/ci) is due to the presence of a real barrier, 
^i]l[o,i]- 

The relation between outgoing wave boundary conditions and wave behaviour at the 
boundary of a bounded absorber is useful for applications and experiments as well as 
for calculations. Experiments to observe acoustic (or electromagnetic) resonances and 
scattering are generally conducted in anechoic chambers, which are lined with baffles of 
sound-absorbing material. These baffles prevent incoming reflected waves from interfering 
with the experiment. Just as one can mimic the "radiation-only" property of an infinite 
domain with a finite absorber, models set in infinite domains are often approximations of 
models over a large finite domain in which the medium through which waves propogate is 
slightly dissipative. 
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